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We estimate the amount of vorticity generated at second order in cosmological perturbation 
theory from the coupling between first order energy density and non-adiabatic pressure, or entropy, 
perturbations. Assuming power law input spectra for the source terms, and working in a radiation 
background, we calculate the wave number dependence of the vorticity power spectrum and its 
amplitude. We show that the vorticity generated by this mechanism is non-negligible on small 
scales, and hence should be taken into consideration in current and future CMB experiments. 
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I. INTRODUCTION 

Cosmological perturbation theory is a crucial tool for the study of the early universe. Linear theory (see, e.g. Refs. [IJ- 
7 ) allows us to model the large scale structure of the universe [5] and calculate the anisotropics of the cosmic microwave 
background (CMB) and compare them to the data from experiments, e.g. Wmap. However, with the advent of new 
data sets, such as those from the full seven year running of Wmap [5] or those anticipated from the recently launched 
Planck satellite |10j . which are superior in both their quality and quantity to previous ones, we are now in a 
position to test predictions of higher order perturbation theory against observational data. Conversely, perturbation 
theory beyond linear order is necessary in order to extract the maximum amount of information from these data sets, 
and thus utilise them to their full potential. The second order theory has been developed by several authors (see 
e.g. Refs. [21 ITTUIrJ] ). and there have been studies of the applications of the theory ranging from calculations of the 
non-Gaussianity of the primordial perturbations to investigations of higher order phenomena produced by couplings 
between lower order perturbations. 

It is well known that at linear order in perturbation theory, different types of perturbations (classified as scalar, 
vector and tensor perturbations, according to their transformation behaviour on three dimensional hypersurfaces) 
decouple [B], However, at second order (see, e.g. Ref. |Hj and references therein), and beyond [T7], this is no longer 
the case, and mode couplings exist which give rise to novel, interesting features. One well-studied example of such a 
higher order effect is that of gravitational waves sourced by first order density perturbations. This topic has received 
a lot of recent attention [TM2T] and the contribution of the second order tensors to the tensor spectrum has been 
studied in some detail. 

The production of vorticity at second order in perturbation theory is another effect that only becomes important 
at this and higher orders in a Friedmann-Robertson- Walker background, and is the focus of the current work. For 
the case of a barotropic fluid Ref. |22| showed that vorticity is not produced at any order in perturbation theory. In 
Ref. [2 3) we showed, dropping this assumption of barotropicity, that vorticity is indeed generated at second order. 
We considered a generic perfect fluid, i.e. a fluid whose energy momentum tensor is diagonal (see, e.g. Ref. [24]). The 
source term of the vorticity evolution equation contains terms coupling the linear energy density and non-adiabatic 
pressure (entropy) perturbations. We derived the evolution equation for the second order vorticity tensor, LJ^ijj 
assuming zero first order vorticity and no anisotropic stress, as 

^ - m& 2ij = {mv 1{i 8P naAltj] + W^IM j , (L1) 

P0+-M) I PO + -U) J 

where a(rf) is the scale factor, H. the Hubble parameter, c 2 the adiabatic sound speed, po an d Pq are the energy density 
and pressure in the background, Vu, 5P aaJ n, $Pi, ar e the covariant velocity, the non-adiabatic pressure, and the energy 
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density perturbation at first order, respectively, and a prime denotes differentiation with respect to conformal time. 
From this equation, we can see that even in the absence of anisotropic stress, terms quadratic in the first order 
perturbations will act as source terms for the vorticity tensor at second order. This is very different from the case 
at linear order where, in the absence of anisotropic stress, the vorticity has no source term and will decay with the 
expansion of the universe |2"5H2"7] . 

To keep our results conservative and our calculation as simple as possible and hence analytically tractable, we 



assume that the source term in Eq. (1.1 1 is dominated by the second term. Then, choosing the radiation era as our 
background in which c 2 = 1/3, the evolution equation simplifies to 

9a 

~ H-Uiii = ^-2 S Pl,U SP nadl,i\ ■ (1-2) 

°Po 

The main goal of this paper is to get a handle on the power spectrum for the second order vorticity sourced by the 



coupling between first order density and non-adiabatic pressure perturbations, using Eq. (1.2 1. We take as input 
power spectra for the source terms simple power laws in time and wavenumber (in Fourier space), and analytically 
solve for the evolution of ui2ij ■ 

The paper is organised as follows. In the next section we introduce our definitions and present the governing 
equations, followed by an analytical solution for the linear energy density perturbation in the presence of entropy. In 



Section III we solve for the evolution of the vorticity. We discuss these results and conclude in Section IV 

In this paper we use conformal time throughout. Latin indices i,j,k, take the value 1,2, or 3, and we denote a 
partial derivative with a subscript comma. The order of the perturbations is denoted with a subscript immediately 
after a perturbed quantity and we work in the uniform curvature gauge. 



II. DEFINITIONS AND EQUATIONS 

The Friedmann-Robertson- Walker (FRW) metric tensor, up to and including second order perturbations has, in 
the uniform curvature gauge, the line element (see, e.g. Ref. [14] ) 

ds 2 = a 2 {r)) [-(1 + 20! + (j) 2 )dr, 2 + (2B U + B 2l )dr)dx l + S^dxHaP] , (2.1) 

where we have neglected tensor perturbations, as they will not change our results qualitatively, and are considering a 
flat (K — 0) background. Here, (j>\ and (f> 2 are the lapse function at first and second order, respectively, and Bu and 
B 2i denote the shear in this gauge. We consider the matter content of the universe to be a perfect fluid, with equation 
of state Pq = wpo, and spatial three velocity Uj. The covariant velocity perturbation, Vi, is defined as Vi = -Bj + i>j 

In the remainder of this section we present the evolution and constraint equations to first order in cosmological 
perturbation theory necessary to derive an analytical solution for the first order density perturbation in the uniform 
curvature gauge. Having found the exact analytical solution we approximate it by a power law valid at early times. 



A. Governing equations 

The evolution and constraint equations are given by energy-momentum conservation and the Einstein field equations, 
respectively. In the background, these are the familiar continuity and Fricdmann equations of a FRW universe, namely 

P ' Q = -m( Po + P ) , (2.2) 

2 8nG 2 ,„ Q x 

ri = —^~ a Po ■ ( 2 - 3 ) 

At first order, the governing evolution and constraint equations are, in the left- and right-hand columns respectively, 

6p' 1 + m(6 Pl + 5P 1 ) + (p + P )V 2 v 1 = 0, Ufa = -AirGa 2 ?! , (2.4) 

V{ + AUVi + (p + P ) <t>i + SPt = , 3H 2 0i + HV 2 B 1 = -^Ga 2 5 P i , (2.5) 

where V 2 = did 1 denotes the spatial Laplacian, and we have introduced the rescaled covariant velocity perturbation, 
Vi, for notational convenience, as 

Vx = (po + -Po) (vi + Br) . (2.6) 
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We refrain from deriving the governing equations Eq. (2.4) and Eq. (2.5) here, and instead point the interested reader 
to, e.g. Ref. [13], for details. Using the Friedmann equation, Eq. (2.3), the right-hand equation in Eq. (2.4) can be 
rewritten as 



3_H 



Vx, 



which can then be used to write the right-hand equation of Eq. (2.5) as 



2 Po 



3H 

2po 



Spi 



Hence the evolution equations then become 

Sp[ + \U (3 + id) (5pi 



ZUSPi - k Vx - \U 2 (1 + w) Vi = , 



Vx + \U (5 - 3w) Vx + SPx = , 



(2.7) 

(2.8) 

(2.9) 
(2.10) 



where we are working in Fourier space, k being a comoving wavenumber (as usual). Equations (2.9) and (2.10) make 
up a system of coupled, linear, ordinary differential equations. Given an equation of state and initial conditions, 
this system can be solved immediately (for a given k), numerically. One can also obtain a qualitative solution by 
considering a system of two equations such as this one. However, if one wants to solve the system quantitatively, 
and analytically, it is ea sier to rewrite the system as a single second order differential equation, which we do in the 
following. We solve Eq. (2.9) for Vx and get 



Vx = v ( Vi +%H(3 + w) 5 Pl + m5P 1 



(2-11) 



where X = k 2 + ^H 2 (1 + w) . After some further algebraic manipulations of Eq. (2.11), and using Eq. (2.10), we 
arrive at the desired evolution equation 

X' 



Sp'l 



7H- 



X 



2 

m5P[ 



H[ w' 
H + {3 + w 
H^ X^ 

n x 



x 1 



i 



H(5- 3w) 



+ ^H 2 (5-3w)+X 



Sp 1 
SPi 



0. 



(2.12) 



Equation (2.12) is a linear differential equation, of second order in (conformal) time. It is valid on all scales and for a 
single fluid with any (time dependent) equation of state. Furthermore, it assumes nothing more than a perfect fluid 
and hence allows for non-zero non-adiabatic pressure, or entropy, perturbations. 



Solutions 



Having derived a general governing equation (2.12) valid in any epoch, we now restrict our analysis to radiation 
domination, where the background equation of state parameter is w = | and the adiabatic sound speed is c 2 = |. 
During the radiation era the scale factor scales as a oc r], and hence H oc 77 -1 . We note that the first order pressure 
perturbation can be expanded as (see, e.g. Ref. [2"8] ~) 

SPx = c 2 6 Pl + (5P nadl , (2.13) 

where the non-adiabatic pressure perturbation is defined, in terms of the energy density, p, and entropy, S, as 



6P 



nadl 



dP 
~dS 



SSx 



(2.14) 



Then, the general governing equation, Eq. (2.12), becomes, in radiation domination, 

m 2 



Sp'l + 4H12 



/-- ( ,^;^ + r 2 + 3 (fc2 + 6 ^ 



72H 4 



-3HSP lna(i ' 



9H 2 + 36 



k 2 + m? 



8px 



k 2 + en 2 



SP, 



nadl 



(2.15) 
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For the case of zero non-adiabatic pressure perturbations the second line in Eq. (2.15) vanishes, and the resulting 
equation can be solved directly using the Frobenius method, to give 



Spi(k,i]) = Ci(fe)r? 



-5 



cos 



h L)kn-2 S in( k X)V3 



V3J 



V3j 



C 2 (k)rf 



2 cos 



krj 

7! 



sin 



kr] \ 



krj 



(2.16) 



where C\ and Ci are functions of the wave-vector, k. For small kr), the trigonometric functions can be expanded in 
power series giving, to leading order, the approximation 



S Pl (k,ri) ~ A{k)ki 1 - i + B(k)r)- 



(2.17) 



for some A(k) and B(k), determined by the initial conditions. 

In order to solve for a non vanishing non-adiabatic pressure, we make the ansatz that the non-adiabatic pressure 
grows as the decaying branch of the density perturbation in Eq. (2.17), i.e., 



(2.18) 



This assumption is well motivated, since we would expect the non-adiabatic pressure to decay faster than the energy 
density. This gives the solution 



8p x (fe, rj) = C\ (fe) T) 5 cos 
+ C 2 (k) v - 5 



2 cos 



■ r> • ( kr] 
kn — 2 sin — = 
V\/3 

Vo + sin 



kr] 

7! 



V3 

krA 
V3j 



kr] 



3D(k)k A r] 



A„-5 



(2.19) 



Anticipating results from the following section, we put the vector dependence of the Sp 1 and <5P na di into Gaussian 
random variables E(k), allowing us to write, for example, 5p 1 (k, rj) = Sp 1 (k, rj)E{k), where k — \k\. Then as a further 
approximation for the scalar fc-dependent source terms, we can expand Eq. (2.19) to lowest order in krj, and use the 
aforementioned ansatz for the non-adiabatic pressure perturbation, Eq. ( |2.18 ). Taking the functions A(k) and D(k) 
to be power laws in k this gives 



where A and D are yet unspecified amplitudes and /3 and a undetermined powers. 



(2.20) 



III. SOLVING THE VORTICITY EVOLUTION EQUATION 



Having obtained a solution for the first order energy density perturbation in the previous section, we can now turn 
towards solving the evolution equation for the second order vorticity. Note that, to avoid notational ambiguities, 
we omit the subscript denoting the order of the perturbation in the following: the vorticity tensor is a second order 
quantity, and the energy density and non-adiabatic pressure perturbation are first order quantities. 



A. The power spectrum of the vorticity 



We now derive the power spectrum for the second order vorticity tensor. Recall from Eq. (1.2) that the vorticity 
evolves, during radiation domination, according to 

9a 

CJy - HuJij = ^-2<fy>,[j£-Piiad,i] = Sij . (3.1) 

The source term can then be written, in Fourier space, as 

S l3 {k,r]) = ~j0j^ J cPkihkj ~ k^SP^ik^Spik - k) , (3.2) 



where we have defined 
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It is easier to work with the vorticity vector, obtained by contracting the vorticity tensor with the totally antisymmetric 
tensor as Wi(x) — Eijkto 3 (x) (see e.g. Ref. [22]). One can define a source vector in an analogous way, which can be 
decomposed in terms of basis vectors as 



Si(x) = 



1 



(2tt)3/2 



A Ak- 



d A kS A {k)e? e 



1 



(2tt) 3 / 2 



d A k 



S\ei 



S 3 h 



Ak-i 



This enables us to write the evolution equation, Eq. (3.1), as 



V) ~ Hu A (k, 7]) = S A (k, rj) . 



(3.4) 



(3.5) 



for each basis state, A. We only have to consider the case A = 1, and therefore drop the subscripts in the following 1 . 
The power spectrum of the vorticity is defined, in the usual manner, as 



(uj*{k 1 ,i])uj(k2,r])) 



2tt 



<S(fci-fc 2 )P u (M) 



(3.6) 



where the asterisk denotes the complex conjugate. Since Eq. (3.5) is a simple first order differential equation in 
conformal time, it can be integrated directly, which then enables us to write the correlator for the vorticity as 



fh 1 



>)0 



ij 2 1 (S* {k 1 ,f]i)S{k 2 ,fi 2 )) dfiidf) 2 . 



(3.7) 



'/o 



where 770 is a constant denoting an initial time. Thus, we can relate the power spectrum, V u {k, 77), to the correlator 
of the source terms by equating Eq. (3.6) and Eq. (3.7). Then, assuming that the fluctuations 5p and <5P na d & re 
Gaussian, we can put the directional dependence into Gaussian random variables E{k) 1 for example, by writing 
Sp(k,r]) — Sp(k,rj)E(k). These Gaussian random variables obey the relationship (E(ki)E(k2)) = <5 3 (fci + k 2 ). Doing 
so enables us to evaluate the correlator of the source terms, (<S*(fci, rf)S(k2, 77)), and obtain, after some algebra, the 
expression for the power spectrum of the vorticity 



"Pu>(k,r)) = 



k 5 rf 
4tt 4 



f(m)f{f)2)Vi % dfjidfk / d^kieikySPn^k^^Spdk + klfn) 



SP na d(\k + k\,fi 2 )5p(k,fi2) - SP nad (k,fi 2 )5p(\k + k\,fj 2 ) 



Substituting the approximations (2.20) into the above gives 

2 



V w (k, rj) = 



rj 



256 4tt 4 



(ADY 



In 



d 3 k{eik l ) 2 k a \k + kf[\k + k\ a k p - k a \k + k\ 



(3.8) 



(3.9) 



where we have performed the temporal integral by noting that as mentioned above, a oc rj during radiation domination, 
and thus po oc 77~ 4 . To perform the /c-space integral, we first move to spherical coordinates oriented with the axis in 
the direction of k. Then, denoting the angle between k and k as 9, the integral can be transformed as 



f d 3 k -> 2vr f " k 2 dk [ 
J Jq Jo 



sm I 



(3.10) 



where the prefactor comes from the fact that the integrand has no dependence on the azimuthal an gle and k c denotes 
a cut-off on small scales. Noting that, in this coordinate system /q<? = fcsin^ the integral in Eq. (3.9) becomes 



I(k) = 2tt 



k c />7T 



^0 



: sin d6 dk sin 2 Qk$ 



l + (fc/fc) 2 + 2(fc/fc)cos6» 



0/2 



(V(l + (k/kf + 2(k/k) cos 9) a ^ 2 k^ - ~k a k p {l + (k/kf + 2(fc//c)cos6»)' 3/2 ) (3.11) 



Note that we make this assumption without loss of generality: of the two other basis states, one is zero from the definition and the other 
gives zero contribution to the power spectrum after the following calculations. 
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Finally, in order to solve this integral we change variables again to dimensionless u and v defined as [TH] (or similarly 

E2I) 



v= u = (k/k) 2 + 2(k/k) cos 6 = y/l + v 2 + 2i 



for which the integral (3.11) becomes 



I(k) = a+ ^+ 5 f kc/k f V+1 uduv 3 dvu^v a {l-^ (u 2 -v 2 -lf) 



J\v-l\ 



(3.12) 



(3.13) 



B. Evaluating the vorticity power spectrum 



In order to perform the integral Eq. (3.131 derived above, we need to specify the exponents for the power spectra 
of the energy density and the non-adiabatic pressure a and j3. The energy density perturbation can be related to the 
curvature perturbation on uniform density hypersurfaces, £, during radiation domination through |14] 



(3.14) 



and hence the initial power spectra can be related as (5pi n iSp- u 
of the initial density perturbation 



16 Po-yini (CiniCini), and we get the power spectrum 



Spini OC Cini OC 



(3.15) 



where ko is the Wmap pivot scale and n s the spectral index of the primordial curvature perturbation [9] . This allows 
us to relate our ansatz for the density perturbation, Eq. (2.201, to the WMAP-data which gives 



Sp = 5 pi, 



k 



1 



— ^.iniPOi] 



k 



(3.16) 



From this, we can read off that P = \ (n s + 1) ~ 1 and the amplitude A — Ai n ipoini- We have some freedom in choosing 
a, however would expect the non-adiabatic pressure to have a blue spectrum, though the calculation demands a ^ j3. 
Using the notation of Ref. [9] we get 



A 2 = pLi^(fco) 2 - pLi*o 6 Ak, 
where we also have the ratio 



D 2 = p 2 ini V s (k Q y = p 2 ini k^A% 



A 2 
A 2 



a(k ) 



and therefore, 



(AD) 2 



1 - a(k ) ' 
a{k 



A 8 4 k -12 



1 - a(k ) 



(3.17) 



(3.18) 



(3.19) 



where we can substitute in numerical values for A 2 ^ and a(ko) from Ref. [5] later on. 
Then, making the choice a = 2, the input power spectra are 



8p(k, 77) = A 



SP nad (k,r]) = D 



for which the integral Eq. (3.131 becomes 



I(k) 



.11 



' u 2 duv 5 dv(l- A;(u 2 -v 2 -I) 2 ) 
|„_i| V Av 2K J 



2 / \ -5 
H 



2 2 
u v — v u 



(3.20) 



(3.21) 
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We can then integrate this analytically, taking care when it comes to evaluating the limits, to give 



I(k) 



16 
135' 



k„ k 



12 
245' 



7;„4 



, k„ k 



5;J3 



1575 



, k„ k 



(3.22) 



which clearly depends upon the small scale cut-off, as expected. For illustrative purposes, let's choose k c = 10 and plot 
the solution I{k) in Figure 1. The left hand figure shows that the amplitude of the integral grows as the wavenumber 
increases. The right hand figure shows a turn around and a decrease in power at some wavenumber (in fact, for a 
non-specific cutoff, this point is at 3.7375fc c ). However, we note that this value is greater than our cutoff, and therefore 
unlikely to be physical. 




- 2x10 11 



1x10 11 




(a)Small range of k < k c . (b)Wide range of k including k > k c . 

FIG. 1: Plots of I(k) for the illustrative choice of k c = 10. 



Then, using the above, and noting that we must use the input to the temporal integrals as 

-4 



>7 



PO — POini 



V 



we obtain the power spectrum for the vorticity, for general k c as 

2 



(3.23) 



•>! A .5 >r 1 ,..o 



r)\ a(fco) 



Pjk - Il] = 256^4^'^ I " ) ( ~ ^ I k ^ A ^ k > 



i]o J \ 1 - a(k ) 



16 k$ f k \ 7 12 k 2 c ( k \ 9 4 (k 
135 14 [koj + 245 fcf V^o / ~ 1575 V^> 



(3.24) 



and on substituting in values from Ref. [9], taking a conservative estimate for a(fco), being 10 per cent of the upper 
bound 



P u (k,ri) = ?7 2 ln 2 f^- 



O.sT lQ- 12 k? ( AVmpc 11 +3.73 x 1Q- I8 fc c 7 ( k 1 



A'o 



Mpc J 



7.71 x lO" 25 /^ f ^- ) Mpc 7 



(3.25) 



and for our above choice of k c — lOMpc 



V„(k,ri) = 77 2 ln 2 ( — 



0.87 x 10" 



_fc_V 
A' 



3.73 x 10 1 I - 7.71 10 



-20 



11 



Mpc 2 



(3.26) 
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IV. DISCUSSION AND CONCLUSION 



In this paper we have obtained the first realistic calculation of the amount of vorticity generated at second order 
in cosmological perturbation theory by allowing for a non-adiabatic pressure perturbation, as would be the case in 
a multi-field or multi-fluid system. By solving the governing equations at first order, we have been able to obtain 
an approximate input power spectrum for the energy density perturbation. Then, making the ansatz that the non- 
adiabatic pressure perturbation has a bluer spectrum than that of the energy density in order to keep the non-adiabatic 
pressure sub-dominant on all scales, we obtain an analytical result for the vorticity. Our results show that the vorticity 
power spectrum has a non-negligible magnitude which depends on the cut off, k c and the chosen parameters. As this is 
a second order effect the magnitude is somewhat surprising. We have also shown that the result has a dependence on 
the wavenumber to the power of at least seven for the choice a = 2. Therefore the amplification due to the large power 
of k is huge, rendering the vorticity not only possibly observable but also important for the general understanding of 
the physical processes taking place in the early universe. 

As discussed above, we have only approximated the input power spectra. This approximation necessarily only holds 
for some time and some range of wavenumbers. The next step will be to solve the system of equations numerically, 
both the vorticity evolution equations and the input power spectra. This can be done for various settings, e.g. by 
calculating <5P na d from multi-field scalar field model [3DJ at or during the end of inflation, or for multi-fluid system 



during later epochs, evolving dp and using Eq. (2.12) 



One prospect for observing early universe vorticity is in the B-mode polarisation of the CMB. Both vector and 
tensor perturbations produce B-mode polarisation, but at linear order such vector modes decay with the expansion 
of the universe. However, vector modes produced by gradients in energy density and entropy perturbations, such as 
those discussed in this work, will source B-mode polarisation at second order. Furthermore, recently it has been noted 
that vector perturbations in fact generate a stronger B-mode polarisation that tensor modes with the same amplitude 
[31 j . Therefore, it is feasible for vorticity to be observed by future surveys. Finally, an important consequence of 
vorticity is the generation of magnetic fields, and so early universe vorticity could play a crucial role in determining 
the origin of primordial magnetic fields. We will investigate this in a future publication [32] . 
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